A kinetic theory of plastic flow in soft glassy materials 
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A kinetic model for the elasto-plastic dynamics of a flowing jammed material is proposed, which 
takes the form of a non-local - Boltzmann-like - kinetic equation for the stress distribution function. 
Coarse-graining this equation yields a non-local constitutive law for the flow, introducing as a key 
dynamic quantity the local rate of plastic events. This quantity, interpreted as a local fluidity, is 
spatially correlated, with a correlation length diverging in the quasi-static limit, i.e. close to yielding. 
We predict finite size effects in the flow behavior, as well as the absence of an intrinsic local flow 
curves. These features are supported by recent experimental and numerical observations. 



Soft amorphous materials such as foams, emulsions, 
granular systems or colloidal suspensions display com- 
plex flow properties at high enough concentrations, in- 
termediate between that of a solid and a liquid: at rest 
they behave like an elastic solid, but are able to flow 
"like a liquid" under sufficient applied stress [U El [3l |4] . 
This mixed fluid/solid behavior occurs above a thresh- 
old volume fraction associated with the appearance of 
a yield stress ad- The yielding behavior makes such 
systems particularly interesting for applications - from 
tooth paste, coatings to cosmetic and food emulsions -, 
but fundamentally difficult to describe 0E1I7]. Further- 
more, it has been recognized over the recent years that 
this yielding behavior is, in most cases, associated with 
peculiar spatial features. This takes the form of inhomo- 
geneous flow patterns, such as shear-bands [TJ |3l [8], 
or cooperativity in the flow or deformation response 
0[TOl[ni[TJ[T3J[Tl[I51[TS], potentially associated with 
non-locality in the constitutive rheological law [TT] and 
dependence of the flow on the nature of the boundaries 
[U [TT] . While such features appear to be generic for this 
class of materials, suggesting a underlying common flow 
scenario, a consistent framework linking the global rhc- 
ology to the local microscopic dynamics is still lacking. 

In this paper, we present a kinetic elasto-plastic (KEP) 
model, which aims at constructing such a link between 
the microscopic and the macroscopic scales. Starting 
from a kinetic elasto-plastic description of the dynam- 
ics, we derive systematically a (non-local) generic consti- 
tutive law for the flow, obtained by coarse-graining the 
microscopic spatio-temporal dynamics. The predictions 
of the KEP model will be shown to capture most features 
of the rheology of yield fluids, and in particular the re- 
cent experimental demonstration of cooperativity in the 
flow behavior of jammed emulsions [TT] . 

The KEP model, which is detailed below, is based on 
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the generic picture which has emerged recently for the 
dynamics of soft glassy materials [T4l ITcH ITU] . In these 
materials, flow occurs through a succession of elastic de- 
formations and localized plastic rearangements associated 
with a microscopic yield stress. These localized events in- 
duce long range elastic modifications of the stress over the 
system, thereby creating long-lived fragile zones where 
flow occurs. Flow in these systems is thus highly cooper- 
ative and spatially heterogeneous: a dynamically active 
region will induce agitation of its neighbours and thus 
a locally higher rate of plastic rearrangements. Correla- 
tions between plastic events are accordingly expected to 
exhibit a complex spatio-temporal pattern [14] . 

The KEP model - Describing these complex dynami- 
cal processes is a formidable task, and to get further in- 
sights, we propose on purpose a schematic model, relying 
on a few simplifying assumptions. To this end, the KEP 
model extends on an approach first proposed by Hebraud 
and Lequeux (HL) [B] , by describing spatial interactions 
between plastic events: the sample is divided into ele- 
mentary blocks i of size a (typically the size of individual 
particles), carrying a scalar shear-stress di [24], and the 
system is described in terms of the block stress distribu- 
tion Pi(a,t). The latter evolves via three mechanisms: 
an elastic response, under an externally imposed shear 
rate 7°; a stress relaxation due to local plastic events; the 
modification of stress due to the plastic events occuring 
in other blocks, transmitted spatially via clastic interac- 
tions. Various simplifying assumptions are made to de- 
scribe these processes. First, the local plastic events will 
be assumed to occur above a a local threshold value of 
the stress a c . The elastic propagation of the shear stress 
is captured using the stress-stress elastic propagator 11^ 
[14"] . relating the stress relaxed at a block i due to the 
occurence of a localized plastic event in another block j: 
b~oi = Ilij • Sa-j where 5a j is the relaxed stress at block 
j (we assume here a full stress relaxation Sa j = —o~j). 
Furthermore, in order to get a closed kinetic equation, 
we propose a decoupling of the plastic-event dynamics, 
in the same spirit as the Boltzmann Stosszahlansatz. Fi- 
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nally, as a first approach, convection is neglected. Al- 
together, within these simplifying hypothesis, the KEP 
equation for Pi(a,t) takes a Boltzmann-like form 

d t Pi(a,t) = -GoWrPfat) - 6(|(7| ~ <7c) Pi((7,t) 

T 

+Ti(t)S{a)+£(P,P) (1) 

with G is the elastic modulus, O the Heaviside function. 
The non-local Boltzmann like operator C(P, P) is defined 
according to : 

[Pjia'^Piia + da^-P^Piia)] (2) 

with 6<7i — Hi j 5a j — —ILija' and 1/r acounts for the 
relaxation rate when the stress is larger than a c . C(P, P) 
describes the gain and loss contributions for the proba- 
bility Pi(a,t) due to events occuring in other blocks, in 
full analogy with the Boltzmann equation. The rate of 
plastic events, Ti(t), entering Eq. 0, is defined as 

r*(*) = / 9(l a '}- ac) P(<r'W, (3) 

In its above form, the KEP equation remains difficult 
to solve analytically. To proceed further, we formally 
expand the Bolztmann operator C(P, P) for small stress 
variations da and retain only the first terms of the ex- 
pansion. The further simplification 8a j ~ —a c (valid 
for small 7) is also made. This simplifies Eq. to a 
Fokker-Planck equation: 

d t Pi(a,t) = -GoWrPfat) - e(M ~ ac) Pi(M) 

r 

+T l (t)S(a) + D l dl 2 P l (a,t) (4) 

In this equation, 73 is the local shear-rate (7^ = 7° + 
I Y^j^i nija c Tj) and the coefficient _D 4 quantifies what 
appears as a stress diffusion induced by the occurence 
of plastic events. A key result is that stress diffu- 
sion is related to the rate of plastic events over the 
whole system via the self-consistency relationship: Di = 
kJ2j^i^ij a c^j' therefore making the Fokker-Planck 
equation non- linear. 

Coming back to continuous spatial variables, a closed 
system of equations is obtained for the local stress diffu- 
sion D(r,t), rate of plastic events r(r,i) and stress dis- 
tribution P(a,r,t). Eq. Q keeps the same form (with 
i — > r), while a small slope approximation of the self- 
constistency equation for D provides a non-local relation- 
ship between stress diffusion and rate of plastic events: 

D(r,t) = mAT(r,t) + aT(r,t) (5) 

with A the spatial Laplacian. In this equation, two key 
parameters have been introduced: a coupling parameter 



a, here defined as a = a% -^f j'i anc ^ an inhomo- 
geneity parameter m = a 2 a 2 H 2 ln , with II n „ the nearest 
neighbour (block-to-block) propagator. In the following 
we will make use of dimensionless variables, t = t/r, 
r = r/a, a = a/a c , 7 loc _= 7i oc G r/a c , m = m/a 2 a 2 , 
a = a/a 2 , f = Tt, and D = Dr/a 2 , but will drop the" 
to simplify notations. 

If inhomogcncitics are neglected [e.g., putting m — 
in Eq. ([EJ], the above set of equations reduce exactly 
to the HL description in Ref. [6]. Anticipating on the 
discussion below, a key result which emerges from the 
HL description is that it predicts a jamming transition 
below a threshold (dimensionless) coupling parameter 
a < a c — I , associated with the building-up of a macro- 
scopic dynamic yield stress, a(j — > 0) = ad- As shown in 
Ref. [B], ad oc (a c — a)^ , with (3 — 1/2 and this dynamic 
yield stress thus quantifies the distance to the jamming 
transition [23 . In the following we shall focus on the 
jammed state, as defined by a non-vanishing ad- 
Constitutive flow rules - In the stationnary state, the 
Fokker-Planck equation, Eq. Q, can be solved ana- 
lytically [B] to give an explicit expression for P(a, r). 
One then deduces the local averaged stress <r(r) = 
/ da' a'P(a' ,r) and the local rate of plastic events T(r) 
(from the normalization condition for P). This provides 
explicit expressions for these quantities in terms of the 
diffusion coefficient D(r) and local shear rate 7(1"). While 
their general expression is rather cumbersome, they sim- 
plify considerably in the quasistatic limit (7 — > 0) and 
close to the jamming point, i.e. small ad- Choosing the 
plastic rate T as the key variable, one obtains the fol- 
lowing expressions in this regime: a = (6T) 1 x 7 and 
D-aT = aia d {a d - a)T + a 2 T 3 / 2 + 0(T 2 ), with a d the 
dynamic yield stress introduced above, and a±, a 2 two nu- 
merical constants [23] . In the following we define / = 6r 
as the fluidity: the latter naturally emerges as intimately 
linked to the rate of plastic events. 

Together with the self-consistency relationship Eq. (|5| , 
relating D to T, these expressions provide a closed set of 
equations. A further linearization allows to rewrite this 
set in the physically meaningful form: 

A/-i(/-/ 6 )=0 (6) 

In this equation we have introduced a "bulk fluidity" 
f b (a): f b (a) = 6(^f ) 2 {a - a d f for a > a d and 
otherwise; and a fluidity correlation length £(ct): £ = 

. — — r for a > a d and £ = . I , m for a < a d . 

These coupled equations constitute the non-local consti- 
tutive flow rules which emerge from the KEP model, and 
are the central result of this work. 

The bulk fluidity fb(a) is the value of the fluidity 
obtained in absence of non local terms, as obtained in 
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FIG. 1: (Top) Dimensionless velocity profiles of a jammed 
material under pure shear flow in a slit of width e = 40 a 
(a = 0.45 < a c ). The full lines are deduced from Eq. Q and 
the dashed lines from the bulk model (m = i.e. £o = 0). 
(a): £o ~ a; from bottom to top a = 1.5,3,4,7 cr^. (b): same 
with £o ~ 8a. (Bottom) Pressure driven flows (c): £o ~ a; 
From bottom to top, the lines correspond to various shear 
stress at the wall a w = 1.5,3,4 ad- (d): same with £o ~ 8a. 
a w — 1.5, 4, 7, 9o~d- The characteristic flow velocity is defined 
as Vb = <7 c e/G T. As a boundary condition, we chose the wall 
fluidity as f w — 2fb(a). 



the HL model: as can be easily verified, it predicts 
a Herschcl-Bulkley expression for the flow rule for low 
shear rates, with ad as the dynamic yield stress: (7(7) = 
ad + Aj n , with n = 1/2 and A a constant depending 
on a. A further key result from Eq. ([6} is the non- 
local nature of the flow curve, which introduces a "flow 
cooperativity length" £. Physically, £ quantifies the spa- 
tial spreading of the plastic activity due to the non-local 
elastic relaxation over the system. Interestingly, the cor- 
relation length diverges at the dynamical yield stress ac- 
cording to £ oc |cr — crrf I 1/ ' 2 oc in agreement with 
recent numerical simulations [12 . 

We emphasize that the non-local flow rule predicted 
by the KEP model, Eqs. is formally identical to the 
cooperative rheology introduced recently to account for 
the flow of confined jammed emulsions |llj . 

Couette and Pressure- driven flows - Let us now dis- 
cuss the solution of this rheological model in various ge- 
ometries. One expects non-local effects to emerge in the 
flow behavior and we introduce a characteristic length 
£0, which we define as £ = £{&d + 5 a), with the some- 
what arbitrary choice 5a — \&d- Non-local effects are 
expected when the size of the system, say e, is compara- 
ble to the correlation length £q- The flow behavior also 
requires boundary conditions for the fluidity at the con- 
fining walls, which - in line with experimental results [TT] 
- we will assume here to be a given function of the stress 
at the walls, f w = f w (a w ). 



For a planar Couette cell made of two parallel walls 
separated by a distance e, the mean shear stress <r(r) is 
spatially homogeneous a = a a and the resolution of Eq. 
^ is straightforward in this geometry. This provides the 
expression of 7: 



7(2) 



fb + (fw - fb) 



cosh[(z- e /2)/£(q )] 
cosh[e/2£(a )] 



(7) 



where z is the distance from the bottom plate. Velocity 
profiles are deduced by integration (assuming here no-slip 
boundary condition at the walls). Figs. [TJa)-(b) show 
the resulting velocity profiles for various characteristic 
lengths £o- In the pressure driven (Poiseuille) geometry, 
the stress varies spatially in the confined channel accord- 
ing to a = VP(z — I ), with VP the constant pressure 
gradient along the slit. The constitutive laws, Eqs. (JfjJ, 
are integrated numerically for various VP and the result- 
ing velocity profiles are displayed in Figs, [ijc)-(d). 

In both Couette and Poiseuille geometries, the flow 
profiles deduced from the non-local constitutive rules de- 
part strongly from the "bulk" prediction (i.e. without 
non-local effects), as soon as the characteristic length 
£0 compares with the confinement [note however that in 
Figs, [ijb)-(d) £ is only £ ~ 0.2 e]. Furthermore the ef- 
fect is more pronounced for the Poiseuille geometry, due 
to the spatial inhomogeneity of the stress map which in- 
deed amplifies the non- locality effect. 

An alternative way of exhibiting non-locality is to plot 
the local flow curve: a(z) versus j(z). Fig. [2] shows the 
result of such a plot for the Poiseuille geometry. As 
evidenced on these curves, the existence of non-locality 
(finite £0) results in a multivalued local flow curve, de- 
parting from the bulk prediction: different values for the 
shear rate 7 are obtained for the same value of the stress, 
obtained here for different pressure gradients VP. A sim- 
ilar multivalued behavior is also obtained upon varying 
the confinement e. In other words cooperativity induces 
finite-size effects in the flow of the jammed material which 
occurs for confinements e ~ £o- Furthermore as seen on 
the local flow curve, Fig. [2] non-local effects do suppress 
the yielding behavior of the fluid. 

Discussion - Altogether these cooperative flow behav- 
ior predicted by the KEP model are in very good agree- 
ment with the recent experimental results for the flow of 
jammed emulsions in microchannels [TT]. In these exper- 
iments, the flow profiles were found to depart from the 
bulk prediction for confinements e typically smaller than 
a few tens of droplet diameters. Furthermore, a consti- 
tutive law similar to Eq. ^ was able to rationalize all 
experimental results, with a cooperativity length scales 
of the order of several droplets diameters. One difference 
however is that no dependence of the cooperativity length 
on shear-rate was reported experimentally, even though 
the flow behavior in the quasi-static limit is difficult to ac- 
cess experimentally and would certainly require a further 
specific investigation. Furthermore similar cooperativity 
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FIG. 2: Local flow curves cr(z) versus 7(2) extracted from 
the velocity profiles in the pressure driven flows of Figs, [ljc)- 
(d), i.e. characterized by two different non-locality parameter 
(a): £0 ~ a; (b): £0 « 8a. The dashed line is the bulk 
flow curve (without non- locality). As expected, deviations 
are more pronounced for the stronger non-locality. Note the 
difference in vertical scales. 



effects are reported in granular flows O [T5j [16] , as well 
as in numerical simulations of deformation of amorphous 
materials [IB] , 

Another prediction of the non-local constitutive law is 
the strong impact of boundary effects on the flow. As can 
be seen e.g. in Eq. |7]), the flow profile within the cell is 
influenced by the wall fluidity. The latter is expected to 
depend on surface properties, e.g. roughness : a smooth 
wall is indeed expected to induce a smaller wall fluidity as 
compared to a rough wall, which in turns will modify the 
shape of the flow profile in the material. The influence 
of boundary roughness on the flow is indeed observed 
experimentally in various systems [H [TTJ [T7] and would 
definitely desserve a more systematic investigation. 

At a more formal level, it is interesting to note that 
the solution of Eq. ^ is the minimum of the square 
gradient 'free-energy': fi(T) = / dr f^(Vr) 2 + u(T, a), 
withw(r,ff) = iaicr d (fj d -CT)r 2 +|a 2 r 5 / 2 + Cl(r 4 ) in the 
limit of small T. This equation is analogous to a Landau 
expansion close to a second order phase transition, with 
the dynamic yield stress as critical point. The rate 
of plastic events T, i.e. the fluidity, plays the role of the 
(dynamic) order parameter. 

Beyond the formal analogy, this suggests a interesting 
alternative point of view for flow inhomogeneities. While 
the present scenario predicts flow inhomogeneities char- 
acterized by a cooperativity length scale, in line with re- 
sults for dense emulsions, a "true" shear-banding would 
merely correspond to a first order phase transition sce- 
nario: i.e. the spatial coexistence between two states 
of different fluidity for the same shear stress. Recent 
experimental findings have connected shear-banding to 
the existence of attractive interactions between particles, 
thereby inducing a flow-structure coupling in the mate- 
rial [20 . The KEP description does not account for these 
features and it would be therefore interesting to include 



local structure variables in the description in order to 
capture such couplings. 

Conclusions - In conclusion, we have derived a non- 
local constitutive equation for the flow of jammed sys- 
tems from a 'microscopic' kinetic elasto-plastic model. 
The resulting description suggests the cooperative nature 
of the flow, in full agreement with recent experimental 
findings [llj . Furthermore, this framework puts forward 
the role of the fluidity as a dynamical order parameter 
characterizing the flow, and here defined as the local rate 
of plastic events in the material. Since one expects plas- 
tic events to trigger local velocity fluctuations, (Sv 2 ), the 
latter quantity could provide an indirect measure of the 
fluidity, in line with granular hydrodynamics approaches 
[2Tj . It is interesting to note that similar observations of 
non-locality have been reported in granular flows close to 
the jamming transition |15[ 116] . suggesting further uni- 
versal characteristics. 
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